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We propose a nonlocal manipulation method to build topological devices through emerging robust 
helical surface states in X 2 = 0 topological systems. Specihcally, in a ribbon of X 2 = 0 Bernevig- 
Hughes-Zhang (BHZ) model with finite-size effect, if magnetic impurities are doped on the top 
(bottom) edge, the edge states on the bottom (top) edge can be altered according to the strengths 
and directions of these magnetic impurities. Consequently, the backscattering between the emerging 
robust helical edge states and gapped normal edge states due to finite-size confinement is also 
changed, which makes the system alternate between a perfect one-channel conductor and a perfect 
insulator. This effect allows us to fabricate topological devices with high on-off ratio. Moreover, it 
can also be generalized to 3D model and more realistic Cd 3 As 2 type Dirac semimetals. 

PACS numbers: 72.25.Dc, 73.20.-r, 73.43.-f, 73.63.-b 


I. INTRODUCTION 

The X 2 = 1 quantum spin hall states and the con¬ 
sequent three-dimensional (3D) strong topological in¬ 
sulators (STI) have attracted great interests since the 
success in their theoretical predictions and experimen¬ 
tal observationsOne of the most important reasons 
is that there exist helical surface states characterized 
by a full insulating gap and protected by the time- 
reversal symmetry^i^ Notably, those surface states are 
spin-momentum locked and exhibit robustness against 
the nonmagnetic impurities, which makes Z 2 = 1 STI 
an ideal candidate for energy-saving topological devices. 

However, there still exist great challenges in the appli¬ 
cation of Z 2 = 1 STI due to their own intrinsic charac¬ 
teristics. Specifically, as a topological device, its on and 
off status should be controlled, which means the system 
can be alternated between a conductor and an insulator 
according to an external field. Since the backscattering 
is suppressed with nonmagnetic impurities, Z 2 = 1 STI 
is born a perfect conductor. Yet for this reason, it is dif¬ 
ficult to realize the off status in Z 2 = 1 STI. By now, 
the most usual manipulation method is to dope mag¬ 
netic impurities or apply an electromagnetic field 
which all encounter the same problems. First, both sur¬ 
faces of the topological system should be doped in or¬ 
der to induce backscattering, which also makes the sys¬ 
tem difficult to return to a conductor. Second, since the 
backscattering only exists where the magnetic impurities 
or external fields are induced, the surface states can only 
be controlled locally. However, experimentally, topologi¬ 
cal insulators usually grow on a substrate, which means 
any manipulation on a substrate surface becomes nearly 
impossible due to the difficulty of doping there. Conse¬ 
quently, a nonlocal method to control the surface states 
close to the substrate in order to build topological devices 
becomes rather an important topic. 

Not long ago, theorists proposed a new kind of Z 2 = 0 


topological system^ii^— Considering finite-size effect, the 
Z 2 = 0 topological system with proper size inherits prop¬ 
erties of Z 2 = 1 STI, besides unique emerging robust 
helical surface states due to the finite-size confinement. 
Here, the words “robust” and “emerging” mean those he¬ 
lical surface states are robust against nonmagnetic dis¬ 
order and only exist in a finite-size system compared to 
the traditional one. Recently, Zhu et al. found the sign 
of these Z 2 = 0 emerging robust helical surface states in 
epitaxial Bi(lll) thin filmsMoreover, in contrast to 
Z 2 = 1 STI, there also exists an unexpected phenomenon 
that the weak anti-localization peak suddenly disappears 
after doping magnetic Co on the top surface of thin Bi 
film. This experiment implies that the surface states on 
the bottom of thin Z 2 = 0 topological system can be ma¬ 
nipulated through a nonlocal method on the top surface, 
which may have great significance in building topological 
devices. 

In this work, we propose a new nonlocal method to 
build topological devices through the emerging robust 
helical edge or surface states in Z 2 = 0 topological 
systems. First, we demonstrate the method in a nar¬ 
row ribbon described by 2D Z 2 = 0 anisotropic BHZ 
model. Without magnetic doping, we find that there ex¬ 
ist one pair of emerging robust helical edge states and 
one pair of gapped normal edge states due to the finite- 
size confinement. Therefore, the system exhibits con¬ 
ducting characteristics. However, if the top (bottom) 
of the 2D BHZ model is doped with magnetic impuri¬ 
ties, the gapped normal edge states on the bottom (top) 
gradually become gapless with increasing magnetization. 
Thus, for moderate magnetization on the top (bottom) 
edge, the system remains good conductor since there is no 
backscattering on the bottom (top) edge. When magneti¬ 
zation becomes strong sufficiently, the Anderson disorder 
can cause strong backscattering on both edges, result¬ 
ing insulating characteristics for the system. Using the 
Landauer-Biittiker formula, we verify the above physical 
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pictures by the two-terminal conductance G, which starts 
from 2e^ jh^ then maintains at the value of e^/h, and fi¬ 
nally falls to zero. This novel phenomenon enables us to 
nonlocally manipulate the on and off status in a Z 2 = 0 
topological system, thus it can also be applied to build 
new topological devices. Our proposal is also extended 
to 3D Z 2 = 0 topological systems, such as Wilson-Dirac 
model. And finally, we further discuss the realization of 
this proposal in Cd 3 As 2 type Dirac semimetals. 

The rest of this paper is organized as follows. In Sec.II, 
we first numerically demonstrate this nonlocal engineer¬ 
ing method in a 2D anisotropic BHZ model and 3D 
anisotropic Wilson-Dirac model. Then, we exhibit two 
detailed plans to build topological devices with Z 2 = 0 
topological systems. Next, in Sec.Ill, we generalize this 
method to more realistic Cd 3 As 2 Dirac semimetal mate¬ 
rials. Finally, a conclusion is presented in Sec.IV. 


II. 2D BHZ MODEL AND 3D WILSON-DIRAC 
MODEL 


We first take the anisotropic BHZ model in square 
lattice as an example^^i2U The four-band tight-binding 
Hamiltonian in the momentum representation reads: 

H{k) 0 
0 H^{-k) 

H{k) = Tz{m — rrix — rriy ^ rrixco^kx ^ rriyco^ky) 

-]-TxVx sin kx + TyVy sin ky (1) 

where m determines the band gap, Vx^y reflects the Fermi 
velocity, and rux^y represents the hopping amplitude be¬ 
tween nearest-neighbor sites along the x, y directions, 
respectively. Tx^y^z are Pauli matrices representing differ¬ 
ent orbits. During our calculation, the values of these pa¬ 
rameters are chosen as below: m = 1.64, = 0.8, = 

1.2^Vx,y = 3. We consider the ribbon geometry, and the 
width along y direction is chosen as Ly = 30a (a is the 
lattice constant) in order to induce finite-size effect. 

In this work, our plan of magnetic doping is to dope a 
soft-magnetic stripe on one edge or surface of the sample. 
Therefore, the strengths and directions of the magnetic 
impurities are entirely determined by external magnetic 
field, and we consider the effect of magnetic impurities in 
X direction by an adding term MxCFx ^ tq on y = 0 edge 
of the ribbon. Actually, even though it will be difficult to 
accurately control the magnetization strength Mx of the 
soft-magnetic stripe in experiment, an external magnetic 
field rotating in x-z direction will also work. Because Mx 
is determined by M cos 0, and the existence of Mz has 
little influence on the final result. 

Without magnetic impurities, the energy band of the 
system is drawn in Fig. 1(a), which is double degenerated 
due to the spin degeneration. There exist two pairs of 
gapless edge states and two pairs of gapped edge states 
locating ai kx = 0 and kx = respectively. When 
the Fermi energy is set Ep = 0.08 as the black dashed 



FIG. 1: (Color online) (a)-(c) stand for calculated energy 
bands in 2D BHZ model with three different strengths of 
magnetic impurities doped in x direction: Mx = 0,1 and 2.1, 
respectively, (d)-(f) represent corresponding probability den¬ 
sity distributions of the edge states for (a)-(c) at Ep — 0.08, 
respectively. Only the left propagating states are plotted here. 


line shows in Fig. 1(a), the corresponding probability den¬ 
sity |T(^)p is drawn in Fig. 1(d). It is obvious that 
the gapless edge states around kx = ^ mostly locate at 
two edges, representing the emerging robust helical edge 
states. Consequently, the edge states around kx = tt rep¬ 
resent the hybridized gapped normal edge states due to 
the finite-size confinement. 

The energy band with magnetic impurity strength 
Mx = 1 is shown in Fig.1(b). We find that one pair 
of the previously double degenerated edge states around 
kx = ^ opens a small gap, while the other one remains 
close. Moreover, one of the originally double degener¬ 
ated energy gaps around kx = tt becomes smaller. In 
Fig. 1(e), we find that the only existing edge states are 
the emerging robust helical ones locating at ^ = 30a. 

Continuing increasing the impurity strength, the above 
tendency exhibits much more explicitly. If the impurity 
strength reaches as strong as Mx = 2.1, the smaller gap 
in Fig. 1(b) around kx = tt becomes nearly gapless in 
Fig. 1(c). In Fig. 1(f), the corresponding probability den¬ 
sity |T(^)p tells us that the only existing edge states are 
those locating around y = 30a. Specifically, the emerging 
robust helical edge state colored by blue mostly locates 
around y = 30a, while the normal edge state colored by 
red extends to inner space, which means it still exhibits 
properties of bulk states. 

In order to quantitatively analysis the physical pic¬ 
tures of energy bands, we further study the transport 
properties of this system, and calculate the conductance 
G of this model with nonmagnetic disorder in a two- 






























3 




FIG. 2: (Color online) (a) is the conductance G vs Mx for 
different Anderson disorder strength w in a. two-terminal sys¬ 
tem. The parameters El — Er = 2.3 for left and right leads, 
Ec — 0.03 for central part, (b)-(d) are the schematic dia¬ 
grams corresponding to the three characteristic stages of the 
conductance in (a) during the process of doping magnetic im¬ 
purities. (e) is the conductance G \s w with fixed Mx- (f) 
is the conductance G vs Mx or M for magnetic impurities 
doped in order or randomly. 


terminal system. The nonmagnetic disorder exists in¬ 
evitably in the realistic system and is modeled by Ander¬ 
son disorder with random potential uniformly distributed 
in [—w/2^w/2]^ where w is the disorder strength. The 
distance between the two leads is 3000a and the Fermi 
energy of the central region is Ec = 0.03. In order 
to simulate metallic leads at the two terminals, their 
widths along y direction and Fermi energies are chosen 
as Nl = Nr = 200a and Er = Er = 2.3, respectively. 

We first show how the conductance varies with the 
strength of magnetic impurities Mx in Fig.2(a). For the 
red line without Anderson disorder, the conductance un¬ 
dergoes a transition from G = 2e^ jh to then back 

to 2e^ jh. The second plateau of G = 2e^/h shakes vig¬ 
orously, because the normal edge states on ^ = 30a still 
have properties of bulk states as shown by the red line 


in Fig. 1(f). The whole transition processes can be un¬ 
derstood with the help of the schematic diagrams from 
Fig.2(b) to 2(d). First, as shown in Fig.2(b), due to the 
finite-size effect, there exist one pair of emerging robust 
helical edge states and one pair of gapped normal edge 
states^ while the gapped normal ones have no contri¬ 
bution to the conductance. Therefore, the conductance 
exhibits G = 2e^ jh at the beginning. Then, in Fig.2(c), if 
we slightly dope a little amount of magnetic impurities in 
X direction on ^ = 0 edge, which will induce the backscat- 
tering between the counter propagating emerging robust 
helical edge states on ^ = 0, only the emerging robust 
helical edge states on ^ have contribution to the 

conductance. Thus, the conductance falls to G = jh 
fast and then remains unchanged. At last, in Fig.2(d), 
when Mx is strong sufficiently, the gapped normal edge 
states on ^ = 0 owns a relatively large gap, and the cou¬ 
pling between the gapped normal edge states on ^ = 0 
and ^ is really weak. Therefore, the originally 

gapped normal edge states on ^ become gapless 

and could make contribution to the conductance. As a 
result, the conductance gets back to G = 2e^/h. 

However, considering the inevitable impurities in real¬ 
ity, the system shows exotic transport properties, which 
is the central result of this paper. If the strength of An¬ 
derson disorder is not zero, as shown by the green, blue 
and brown lines in Fig.2(a), the conductance never re¬ 
turns to G = 2e^//i, but continues falling to nearly zero 
with strong Mx- This difference can be explained by the 
backscattering between the emerging robust helical edge 
states and the normal ones on ^ = 30a caused by An¬ 
derson disorder. The existence of this backscattering is 
guaranteed by the physical meaning of Z 2 topological in¬ 
variant: gapless helical edge states with even number are 
not robust against disorder. Besides, for the green line, 
when the impurity strength approaches Mx = 1.5, we 
find the conductance exhibits a minimum value during 
the process of decreasing. This critical point is actually 
also the transition point from G = jh to 2e^/h with¬ 
out disorder, which means the Fermi energy just touches 
the bottom of the closing energy band on = tt. In 
this case, the backscattering in gapped kx = channel 
is permitted. Therefore, there exist more backscatter¬ 
ing channels shown in Fig.2(d), and the backscattering 
between the emerging robust helical edge states and the 
gapped normal ones on ^ = 30a is somehow enhanced. 
As a result, the conductance at this point exhibits smaller 
value than the surrounding points. 

Next, let’s fix the magnetic impurity strength M^,, and 
study how the the conductance G varies with Anderson 
disorder strength w. In Fig.2(e), it is obvious that the 
conductance with Mx = 0.5 exhibits stronger robustness 
than that with Mx = 2, which is consistent with Fig.2(a). 
Moreover, for instance, if the disorder strength is chosen 
as re = 2, the “on” status of the red line with Mx = 0.5 
exhibits G « 0.S7e‘^ jh, while the “off” status of the blue 
line with Mx = 2 exhibits G « 0.03e^//i. Therefore, the 
on-off ratio can be achieved about 30 : 1 in this case. 
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Till now, we find that by altering the strength of mag¬ 
netic impurities on one edge, a Z 2 = 0 BHZ model with 
emerging robust helical edge states can be tuned from 
a perfect one-channel conductor to a perfect insulator, 
which can be applied to build topological devices. 

However, in some cases, it may not be convenient to 
alter the external field, and the strength of magnetic 
impurities is fixed. Therefore, let’s consider a sit¬ 
uation that the direction of doped magnetic impurities 
is not ordered but in a random arrangement. For sim¬ 
plicity, we just assume that M is either in x or z di¬ 
rection. During our calculation, 

^ [ ^totah ^totai]-) Q-iid Mfotai ^ The red 

line of Fig.2(f) shows how the conductance varies with 
magnetic impurity strength M when re = 1.5. To have a 
better comparison, a similar re = 1.5 situation with mag¬ 
netic impurities ordered in x direction is also plotted by 
a blue line. As we can see, when M = 2, the conductance 
of the red line still maintains at G ~ O.Se^/h compared 
with G ~ O.le^/h of the blue line when = 2, which 
can be regarded as the on and off status of the current. In 
experiment, this kind of random can be realized in many 
ways, such as by heating and so on. Therefore, topologi¬ 
cal devices built by Z 2 = 0 2D BHZ model with emerging 
robust helical edge states can also be tuned from a perfect 
one-channel conductor to a perfect insulator by altering 
the directions of magnetic impurities on one edge. 



FIG. 3: (Color online) (a): The conductance G vs Mx for 
^2 = 1 STI with m = 1 in BHZ model, (b): The conductance 
G vs Mx for Z 2 = 0 topological system with m = 1.64 in 
BHZ model, (c): The conductance G vs Mx for Ly — 200a 
without finite-size effect in BHZ model, (d): The conduc¬ 
tance G vs disorder strength w in BHZ model for Ly = 200a 
without finite-size effect. In (c) and (d), all parameters of the 
Hamiltonian are chosen the same as Fig.2(a) except Ly. 

For further showing the advantages of the above pro¬ 
posals, we also consider a 2D Z 2 = 1 BHZ model by 
changing the parameter m in Eq(l) from 1.64 to 1 and 
Ly still maintains 30a. As shown in Fig.3(a), no mat¬ 
ter whether there exists Anderson disorder or not, the 
conductance never falls to zero, but keeps G = e^/h all 


the time. Compared with the Z 2 = 0 topological system 
in Fig.3(b), the Z 2 = 1 STI can not turn off the cur¬ 
rent, which means it is not suitable to build topological 
devices. Moreover, we also calculate a Z 2 = 0 2D BHZ 
model without finite-size effect by expanding Ly from 30a 
to 200a. The red line in Fig.3(c) represents how the con¬ 
ductance G varies with magnetic impurity strength Mx 
without Anderson disorder, and there exists two charac¬ 
teristic quantum plateaus G = jh and G = 2e^/h. 
In Fig.3(d), we show how the conductance on these two 
plateaus varies with Anderson disorder strength w. The 
red line and the blue line represent Mx = 0.1 and Mx = 1, 
respectively. Unlike Fig.2(e), the main feature of Fig.3(d) 
is that the conductances of these two lines start decreas¬ 
ing at the beginning, and fall to zero nearly with the same 
Mx. This feature tells us that one cannot get high on-off 
ratio for any disorder strength w or magnetic impurity 
strength Mx- The reason is that there are always two 
pairs of gapless helical edge states around y = Ly, and 
the Anderson disorder can cause the backscattering be¬ 
tween them (the helical edge states around kx = 0 are no 
longer emerging robust ones now). Therefore, if we want 
to switch on or off the current in order to build topologi¬ 
cal devices just by tuning the strength of doped magnetic 
impurities on one edge, the existence of Z 2 = 0 emerging 
robust helical edge states due to finite-size confinement 
is the necessary condition. 

Similar phenomenon also exists in 3D Z 2 = 0 topolog¬ 
ical systems, and we take the anisotropic Wilson-Dirac 
type model as an example The Hamiltonian in a 
cubic lattice reads: 

^{{k) = m{k)ao 0 Tz Va sin kaCTa 0 Tx 

Oi 

m(k) = m + E ma{coska-l) (2) 

Oi 

where a are Pauli matrices in spin space, and a = x,y, z. 
Parameters r,m,Va,ma have the same meanings as the 
2D model of Eq(l). The width along z direction Lz is cho¬ 
sen as Lz = 15a, and the magnetic impurities are doped 
in z direction on 2 : = 0 side. During the calculation, re¬ 
lated parameters are chosen as below: m = 2.26, rrix = 
6.9,my = l.l,mz = 6.S,Vx,y,z = 1-5. In Eig.4(a)-(c), we 
show how the 3D energy band varies with the strength of 
doped magnetic impurities Mz = 0,1 and 2, respectively. 
There exist two nonequivalent Dirac points A(G) and 
B{D), and the blue and pink curves represent the two 
energy bands nearest to the Dirac points. In Eig.4(a), 
the energy bands locating around A{G) are gapless, be¬ 
cause they are the emerging robust helical surface states. 
However, the appearance of the gaps on B{D) tells us 
these surface states are the gapped normal ones due to 
finite-size confinement. In Eig.4(b), Mz = 1, we find that 
the pink curves on A{C), which were gapless at first, now 
open a small gap, because there exists the backscattering 
between the counter propagating emerging robust helical 
surface states on 2 : = 0 induced by the magnetic impu¬ 
rities. At last, in Eig.4(c), when the strength of mag- 
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FIG. 4: (Color online) (a)-(c) correspond to the surface energy bands of Wilson-Dirac model with different strengths of magnetic 
impurities: Mz = 0,1 and 2, respectively. As Mz increases, one of the two double degenerated gapless surface states on A(C) 
opens a small gap, while one of the two double degenerated gapped surface states on B{D) gradually becomes gapless. 


netic impurities is strong sufficiently, all blue curves at 
A{C) and B{D) become gapless. Briefly, Dirac points 
A{C) and B{D) correspond to kx = 0 and kx = tt in 2D 
BHZ model, respectively. And all the physical pictures 
of Fig.4(a)-(c) are exactly the same as those shown by 
Fig.l(a)-(c), respectively. 

By now, we have afforded two kinds of nonlocal meth¬ 
ods to build topological devices in Z 2 = 0 topological 
systems. The first one is to alter the strengths of the 
doped magnetic impurities. Specifically, we first dope 
a soft-magnetic film at the top surface of a thin ^2 = 0 
topological system grown on a substrate. Then, if we only 
apply a weak magnetic field on the soft-magnetic film, the 
whole device is a good conductor. However, if the soft- 
magnetic film is magnetized strong sufficiently, the whole 
device turns into an insulator. The second method is to 
alter the directions of doped magnetic impurities. If the 
doped magnetic film is magnetized in order with moder¬ 
ate strength, the whole device is an insulator. However, 
if the magnetic direction of doped magnetic film is in a 
random arrangement, such as by heating and so on, the 
whole system gets back to a conductor. Since the above 
manipulation methods are all nonlocal ones on one edge 
or surface of the sample, it exhibits great advantages to 
build topological devices through emerging robust heli¬ 
cal surface states in Z 2 = 0 topological systems compared 
with traditional Z 2 = 1 STL 


III. EXTENSION TO REAL MATERIALS 


perimental realization of Dirac semimetal Cd 3 As 2 ^~'— 
we consider a Cd 3 As 2 stripe whose effective low energy 
Hamiltonian H{k) reads as)2S 


H{k) = eo{k) + 


/ M{k) 
Ak- 
Dk+ 
B{k) 


Ak+ 

-M{k) 

B{k) 

0 


Dk- 

B*{k) 

M{k) 

-Ak+ 


B*{k) 

0 


-Ak- 

-M{k) 


(3) 


where eo(fe) = Co + Cifc^ + C 2 {kl + ky), k± = k^ ± iky, 
and M{k) = Mq — Mik^ — M 2 {k^ + ky) with parameters 
Mo, Ml, M 2 < 0. Since B{k) = {akz + f3D)k‘^, B{k) can 
be neglected compared with Dk± term if we only consider 
the expansion up to 0{k^)^^ 



Based on above two simple models, we have in prin¬ 
ciple shown the nonlocal manipulation methods to build 
topological devices with emerging robust helical surface 
states in Z 2 = 0 topological systems. In this section, 
let’s further extend our proposal to more realistic ma¬ 
terials. Though there are several 3D realistic material 
candidateswe restrict our discussion on 2D cases to 
avoid the huge computational requirements. 

Motivated by the recent theoretical prediction and ex- 


FIG. 5: (Golor online) (a)-(c) stand for calculated energy 
bands in Gd 3 As 2 for three different strengths of magnetic im¬ 
purities doped in x direction: Mx = 0, 0.2 and 0.9, respec¬ 
tively. (d)-(f) are the probability density distributions which 
has been integrated along z axis for (a)-(c) at Ep = 0.014, 
respectively. In (d)-(f), only the left propagating states are 
plotted. 

In this paper, Cq^Ci and C 2 are chosen zero, which 
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can simplify our numerical calculation but has no influ¬ 
ence on the physical properties of this model. In fact, 
main results got from this approximation satisfy realistic 
Cd 3 As 2 materials, though they are not quite true in ex¬ 
periment. Other related parameters are chosen as below: 
Mo = -0.44, Ml, 2 = -0.5, A = 1.2, D = 0.55. Consider¬ 
ing 3D stripe geometry in which Lz = 6a and Ly = 20a, 
the above k • p Hamiltonian can be rewritten in tight- 
bind form, and the energy bands of Cd 3 As 2 are drawn in 
Fig.5(a)-(c). Compared with Fig.l(a)-(c), the two Dirac 
points kx = 0 and kx = tt in BHZ model both locate at 
kx = 0 with two inverted bandsIn Fig.5(a), there is 
no doped magnetic impurity. The two pairs of gapless 
surface states represent the emerging robust helical ones, 
and the other two with a visible gap represent the gapped 
normal ones. Taking the Fermi energy ns Ef = 0.014, we 
obtain the distribution of the emerging robust helical sur¬ 
face states in Fig.5(d), and it has been integrated along 
2 : axis. As we can see, the red line and blue line seem 
degenerated, which originates from the coupling between 
the emerging robust helical surface states on ^ = 0 and 
y = 20a. Next, a little amount of magnetic impurities in 
X direction {MxCFx C)to) are doped on ^ = 0. In Fig.5(b), 
we find that one of the two pairs of emerging robust heli¬ 
cal surface states opens a small gap, while one of the two 
pairs of gapped normal surface states tends to close. In 
this case, the only existing gapless surface states are the 
emerging robust helical ones concentrating on ^ = 20a 
as shown in Fig.5(e). At last, when the strength of mag¬ 
netic impurities is strong enough, such as Fig.5(c) with 
Mx = 0.9, there exist two pairs of gapless surface states 
consisting of one pair of emerging robust helical surface 
states and one pair of normal surface states. The proba¬ 
bility density distribution in this case is shown in Fig.5(f). 
The blue and red lines represent the emerging robust he¬ 
lical surface states and the normal ones, respectively. We 
find only the surface states concentrating around y = 20a 
exist now. In fact, every characteristic in Fig.5(a)-(f) 
of Cd 3 As 2 owns its correspondence in Fig.l(a)-(f) of 2D 
BHZ model, and the physical pictures between these two 
models are exactly the same. 



FIG. 6: (Color online) (a): The red and blue lines stand 
for the conductance G vs Mx in Cd 3 As 2 for two Anderson 
disorder strengths rc = 0 and w = 0.5, respectively, (b): The 
red and blue lines stand for the conductance G vs rc in Cd 3 As 2 
for two fixed magnetic impurity strengths Mx =0.2 and 0.8, 
respectively. 


Finally, we also calculate the conductance of Cd 3 As 2 
with Anderson-type disorder in a two-terminal system. 
The Fermi energies of the two leads and the central part 
are El = Er = 0.5 and Ec = 0.01, respectively. With¬ 
out loss of generality, the disorder is only added on the 
y = 20a boundary^^ Fig.6(a) tells us how the conduc¬ 
tance varies with the strength of magnetic impurities Mx. 
In the red line, there is no Anderson disorder and the 
distance between the two leads is 4 x 10^a. Similar to 
2D BHZ model, the conductance falls to a plateau of 
G = e^/h from the beginning G = 2e^/h rapidly, and 
then goes back to G = 2e^/h. However, if the disorder 
strength is chosen w = O.b shown by the blue line, the 
conductance decreases from G = 2e^/h to G = e^/h, 
and never goes back to G = 2e^/h but continues falling 
to nearly zero. Similar as the green line in Fig.2(a), a 
critical point around Mx = 0.6 also exists, but exhibits 
much more explicit than that in BHZ model. Accord¬ 
ing to the conductance behavior under moderate disorder 
strength w = 0.5, we can claim the two-terminal device 
is in “on” status when Mx G [0.1, 0.3] and in “off” sta¬ 
tus when Mx >0.55. In order to better understand the 
performance of the topological device, we also show how 
the conductance varies with Anderson disorder strength 
w with fixed magnetic impurity strength Mx in Fig.6(b). 
When the device is in the claimed “off” status (such as 
Mx = 0.8), the conductance G quickly falls to zero with 
the increasing disorder strength w. In contrast, when the 
device is in the claimed “on” status (such as Mx = 0.2), 
the conductance nearly maintains at G = e^/h, no mat¬ 
ter how the disorder strength w or the length of the strip 
changes. According to Fig.6(b), a relatively high on-off 
ratio about 10 : 1 is obtained when w > 0.4, which is suit¬ 
able for the topological devices. Moreover, the stronger 
the disorder strength is and the longer the stripe is, the 
higher the on-off ratio can be obtained. Nevertheless, 
disorder strength can not be too strong in case of the 
broken of the “on” status. In fact. Fig.6(b) is consistent 
with Fig.2(e), which is also the result of finite-size effect 
and the key to switch on or off the current by doping 
magnetic impurities on one surface. 

To summarize, all conclusions obtained from 2D BHZ 
model can also be found in Cd 3 As 2 type Dirac semimet¬ 
als, which generalizes our proposals building new topo¬ 
logical devices to a much broader prospect. 


IV. CONCLUSION 

In conclusion, we show that in Z 2 = 0 topological sys¬ 
tems, one can switch on and off the current nonlocally 
just by doping magnetic impurities on one edge or sur¬ 
face of the sample, which enables us to build topological 
devices through emerging robust helical surface states in 
Z 2 = 0 topological systems. This proposal is first demon¬ 
strated in 2D Z 2 = 0 anisotropic BHZ model and 3D 
Wilson-Dirac type model in detail. Notably, it is also 
generalized to realistic Cd 3 As 2 type Dirac semimetals. 
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Since the manipulation methods proposed here are all 
nonlocal ones, Z 2 = 0 topological system with emerg¬ 
ing robust helical surface states exhibits much brighter 
prospects in building topological devices than the tradi¬ 
tional Z 2 = 1 STL 
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